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Based on the invariance principle of differential equations, a 
simple adaptive-feedback scheme is proposed to strictly synchronize 
almost all chaotic systems. Unlike the usual linear feedback, the 
variable feedback strength is automatically adapted to completely 
synchronize two almost arbitrary identical chaotic systems, so 
this scheme is analytical, and simple to implement in practice. 
Moreover it is quite robust against the effect of noise. The famous 
Lorenz and Rossler hyperchaos systems are used as illustrative 
examples. 
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Since it was shown in [1] that for some chaotic systems the synchronization is pos- 
sible, synchronization of (unidirectionally) coupled chaotic systems and its potential 
applications in engineering have been a field of great interest over a decade, see 
[2-4] and references cited therein. Due to the different applications, various specific 
synchronization schemes have been proposed in the literature, see [3] and references 
cited therein. However, just as what stated in [5], despite the large amount of effort 
many key issues remain open. One of the central questions is: Given two arbitrary 
identical chaotic systems, how can one design a physically available coupling scheme 
that is strictly guaranteed to produce stable identical synchronization motion (i.e., 
high-quality synchronization)? In most of the rigorous results based on the Lya- 
punov stability or the linear stability, the proposed scheme is very specific, but also 
the added controller is sometimes too big to be physically practical. One practical 
scheme is the linear feedback. However, in such technique it is very difficult to find 
the suitable feedback constant, and thus numerical calculation has to be used, e.g., 
the calculation of the conditional Lyapunov exponents. Due to numerical calcula- 
tion, such scheme is not regular since it can be applied only to particular models. 
More unfortunately, it has been reported that the negativity of the conditional Lya- 
punov exponents is not a sufficient condition for complete chaotic synchronization, 
see [6]. Therefore the synchronization based on these numerical schemes can not be 
strict (i.e., high-qualitative), and is generally not robust against the effect of noise. 
Especially, in these schemes a very weak noise or a small parameter mismatch can 
trigger the desynchronization problem due to a sequence of bifurcations [7] . 

Actually, this open problem, although significant for complete chaos synchroniza- 



tion, is very difficult and cannot admit the optimization solution [3]. For example, 
in [5] a rigorous criteria is presented to guarantee linearly stable synchronization 
motion, but the criteria is so complicated that the specific numerical calculation 
is necessary for particular examples in practice. The similar problem was ever ad- 
dressed in [8]. 

In this letter, we give a novel answer to the above open problem. We prove 
rigorously by using the invariance principle of differential equations [9] that a simple 
feedback coupling with the updated feedback strength, i.e., an adaptive-feedback 
scheme, can strictly synchronize two almost arbitrary identical chaotic systems. 

Let a chaotic (drive) system be given as 

x = f(x), (1) 

where x = (x u x 2 , • • • , x n ) G FT, f(x) = (f 1 (x)J 2 (x) r --J n (x)) : FT -> R n is a 
nonlinear vector function. And let H C R" be a chaotic bounded set of (1) which is 
globally attractive. For the vector function f(x), we give a general assumption. 

For any x = (x 1 , x 2 , ■ ■ ■ , x n ), x = x 2 , ■ ■ ■ , x°) G Q, there exists a constant 
I > satisfying 

I fi(x) ~ fi(x ) \< Imaxj | Xj - x° |, i = 1, 2, • • • , n. (2) 

We call the above condition as the uniform Lipschitz condition, and / refers to 
the uniform Lipschitz constant. Note this condition is very loose, for example, the 
condition (2) holds as long as §^(i, j — 1, 2, • • • , n) are bounded. Therefore the class 
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of systems in the form of (l)-(2) include all well-known chaotic and hyperchaotic 
systems. Consider the variables of (1) as coupling signals, the receiver system with 
variables y G R™ is given by the following equations: 

y = f(y) + e(y-x), (3) 

where the feedback coupling e(y — x) = (eiei, £262, • • • , e n e n ), e$ = (yi — Xi),i = 
1, 2, • • • , n denoting the synchronization error of (1) and (3). Instead of the usual 
linear feedback, the feedback strength e = (ei, e 2 , • • • , e„) here will be duly adapted 
according to the following update law: 

Q = -7ie-,i = l,2,---,n, (4) 

where 73 > 0, i — 1, 2, • • • , n, are arbitrary constants. For the system of 2n equations 
(which is formally called as the augment system for the convenience below) , consist- 
ing of the error equation between (1) and (3), and the equation (4), we introduce 
the following nonnegative function 

V=l±e> + l±he i + L)>, (5) 

z i=l z i=l li 

where L is a constant bigger than nl, i.e., L > nl. By differentiating the function V 
along the trajectories of the augment system, we obtain 

n n 1 

= E eiiMv) - h{x) + e t e z ) - + L)e\ < (nl - L) ]T e\ < 0. (6) 

i=l i=l i=l 

where we have assumed x,y E £1 (without loss of the generality as f2 is globally 
attractive ), and used the uniform Lipschitz condition (2). It is obvious that V = 



if and only if = 0, i — 1, 2, • • • , n, namely the set E — {(e, e) e R 2n : e = 0, e = 
e G R n } is the largest invariant set contained in V = for the augment system. 
Then according to the well-known invariance principle of differential equations [9], 
starting with arbitrary initial values of the augment system, the orbit converges 
asymptotically to the set E, i.e., y — > x and e — > e as t — > oo. 

Obviously, such identical synchronization motion is strict (i.e., high-qualitative), 
global (as long as the chaotic attractor is globally attractive), and nonlinear stable. 
In particular the nonlinear global stability implies that such chaos synchronization 
is quite robust against the effect of noise, namely under the case of presenting a 
small noise the synchronization error eventually approaches zero and ultimately 
fluctuates around zero wheresoever the initial values start. In addition, we note 
that in order to reach the synchronization the variable feedback strength e will be 
automatically adapted to a suitable strength e depending on the initial values, 
which is significantly different from the usual linear feedback. As is well known, 
in the usual linear feedback scheme a fixed strength is used wheresoever the initial 
values start, thus the strength must be maximal, which means a kind of waste in 
practice. However the final strength in the present scheme depends on the initial 
error, thus the strength is generally equal or less than those used in the constant 
gain schemes. Moreover, the present control scheme does not require to determine 
numerically any additive parameters, and is simple to implement in practice since 
the technique is similar to the well-known self-adaptive controller in the control 
theory. Note that although the converged strength is not bigger than the linear one, 
theoretically the strength may be too big to be practical. However the flexibility 



of the strength in the present scheme can overcome this limitation once such case 
arises. For example, suppose that the feedback strength is restricted not to exceed 
a critical value, say k. In the present control procedure, once the variable strength 
e exceeds k at time t — t , we may choose the values of variables at this time as 
initial values and repeat the same control by resetting the initial strength e(0) = 0. 
Namely one may achieve the synchronization within the restricted feedback strength 
due to the global stability of the present scheme, which is slightly similar to idea 
of OGY control [10]. This excellence is absent from the usual feedback scheme. In 
addition, in the present scheme the small converged strength may be obtained by 
decreasing suitably the value of parameter 7, which governs the rate of increasing 
of the feedback strength. As for the question on how the parameter 7 affects the 
convergence rate (i.e., transient time) and the final coupling strength remains to be 
further investigated. 

Note that in the proposed scheme the coupling all the variables may be redundant 
to achieve synchronization, because we find from my proof that it is not necessary, 
e.g., one may set q = (i.e., cancelling the corresponding coupling) if | e$ |<| 
ej |. Actually this case exists in general due to the nonhyperbolicity of chaotic 
attractor, namely near a nonhyperbolic point the divergence rate of trajectories 
(resp. the contraction and stretching of phase space along trajectories) in some 
directions vanishes (or is very small), see the following examples. Although we can 
not give a rigorous criteria to determine only which variables can be used as the 
coupling signals for the general examples, in a concrete model one may determine it 
by the numerical technique, e.g., the calculation of the Lyapunov exponent respective 
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to each direction. Of course for the low dimensional systems the optimal coupling 
variables may be found by directly testing again and again. Especially, if based 
on the calculation of the conditional Lyapunov exponents a chaotic system can be 
synchronized by linearly coupling a variable, then the coupling this variable in the 
present scheme can surely achieve synchronization, see the first example below. 

Next we will give two illustrative examples. Consider the Lorenz system, 

Xi = (3(x 2 — Xi), x 2 = olx\ — xix 3 — x 2 , x 3 = X\X 2 — bx 3 . (7) 

The corresponding receiver system is: 

yi = /%2-J/i) + eiei, 2/2 = ayi ~ Vm - 2/2 + e 2 e 2 , y 3 = ym - by 3 + e 3 e 3 (8) 

with the update law (4). Now let (3 — 10, a — 28, b — |. It has been well known that 
a suitable linear feedback in the second component may synchronize the Lorenz sys- 
tem, so we let €\ = e 3 = (i.e., the time series of only the variable x 2 are selected as 
the driven signal), and set 72 = 0.1. The corresponding numerical results are shown 
in Figure 1 and Figure 2, where the initial feedback strength is set as zero. Figure 1 
shows the temporal evolution of synchronization error between (7) and (8) , and the 
variable strength e 2 . Figure 2 shows that when an additive uniformly distributed 
random noise in the rang [—1,1] (i.e., a noise of the strength 1) is present in the 
signal output x 2 of (7), the synchronization error eventually approaches zero and 
ultimately fluctuates slightly around zero, meanwhile the variable feedback strength 
is affected slightly and can not stabilize. 

As the second example, we consider the Rossler hyperchaos system: 

xi = —x 2 - x 3 , x 2 = xi + 0.25x 2 + £4, 



X3 = 3 + xiX3, Xi = —0.5x 3 + 0.05x4. (9) 
The receiver system is: 

Vi = ~V2 ~ 2/3 + eiei, 2/2 = 2/1 + 0.25y 2 + 2/4 + e 2 e 2 , 

2/3 = 3 + yi2/3 + e 3 e 3 , 2/4 = -0.5y 3 + 0.05t/ 4 + e 4 e 4 (10) 

with the update law (4). Let e± = (we speculate that this is not optimal), 7, = 
l,z = 2,3,4, and initial feedback strength be (0,0,0), numerical results for two 
different cases are shown in Figure 3 and Figure 4, respectively. Figure 3 shows 
the hyperchaotic synchronization and Figure 4 the slight effect of a noise with the 
strength 0.1 which is simultaneously added to the signals x 2 , £3 and X4. 

The above numerical examples show that chaotic or hyperchaotic synchroniza- 
tion can be quickly achieved by the present controller (i.e., the transient time to 
synchronization is very short). In addition, we find from these examples that such 
synchronization is robust against the effect of noise, namely if the expected synchro- 
nization means that the synchronization error is eventually smaller than a threshold 
value (not tends to zero), then the present scheme is physically feasible in the noisy 
case. Moreover, by comparing the converged feedback strength and the correspond- 
ing feedback signals we find that the coupling is indeed small in the examples. In 
addition by testing the other chaotic systems including the Rossler system, Chua's 
circuit, and the Sprott's collection of the simplest chaotic flows we find that the 
coupling only one variable is sufficient to achieve identical synchronization of a 
three-dimensional system. 

In conclusion, we have given a novel answer to the open problem in the field 
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of identical chaos synchronization. In comparison with the previous methods, the 
proposed scheme supplies a simple, analytical, and (systematic) uniform controller 
to synchronize strictly two arbitrary identical chaotic systems satisfying a very loose 
condition. The technique is simple to implement in practice, and quite robust against 
the effect of noise. The control idea may be also generalized to the case of the discrete 
chaotic systems. We also believe that such simple synchronization controller will be 
very beneficial for the applications of chaos synchronization. Especially the similar 
control scheme has been successfully used to stabilize the chaotic neuron model 
[11], and hence the proposed adaptive-feedback synchronization controller can be 
used to explore more reasonably the interesting dynamics found in neurobiological 
systems, i.e., the onset of regular bursts in a group of irregularly bursting neurons 
with different individual properties [12]. 
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FIG.l. The chaos synchronization between (7) and 



(8) is achieved by only the signal X2, where (a)-(c) 



show temporal evolution of the synchronization er- 



ror and (d) shows the evolution of the correspond- 



ing feedback strength e 2 - Here the initial values of 



(x, y, e) is set as (2, 3, 7, 3, 4, 8, 0). 



1 1 




5 10 15 20 25 30 35 40 45 50 

10 1 1 1 1 1 1 1 1 1 1 1 




5 10 15 20 25 30 35 40 45 50 




5 10 15 20 25 30 35 40 45 50 



timet 

FIG. 2. (a)-(d) show the effect on the synchroniza- 
tion error and the feedback strength in Figure 1 when 
a noise of the strength 1 is present in the signal out- 
put x 2 of (7). 
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FIG. 3. (a)-(g) show the hyperchaotic synchroniza- 



tion between systems (9) and (10), and temporal evo- 



lution of the corresponding feedback strength, where 



only three variables of (9), x i} i — 2,3,4 are se- 



lected as the driven signals. Here the initial values 



of (x, y, e) is set as (2, 3, 7, 10, 3, 4, 5, 11, 0, 0, 0) 



1 o 



a)™ 



10 15 20 25 




_25 30 

^ (c) 



_J5 20 



_25 30 



10 15 20 25 



Ji 20 25 



-2 
1 

-2 



J5 20 



10 15 20 25 

time 1 



FIG. 4. (a)-(g) show the effect on the synchroniza- 



tion error and the feedback strength in Figure 3 when 



a noise of the strength 0.1 is simultaneously added 



to the signal outputs x 2 , x 3 and x A of (9). 
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